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Abstract 
Admixture between species is a cause for concern in wildlife management. Canids 
are particularly vulnerable to interspecific hybridisation, and genetic admixture has 


shaped their evolutionary history. Microsatellite DNA testing, relying on a small 
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number of genetic markers and geographically restricted reference populations, has 


identified extensive domestic dog admixture in Australian dingoes and driven con- 
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servation management policy. But there exists a concern that geographic variation 
in dingo genotypes could confound ancestry analyses that use a small number of ge- 
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and then carry out comparisons to domestic dogs. We then perform ancestry model- 
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ling and biogeographic analyses to characterise population structure in dingoes and 
investigate the extent of admixture between dingoes and dogs in different regions of 
the continent. We show that there are at least five distinct dingo populations across 
Australia. We observed limited evidence of dog admixture in wild dingoes. Our work 
challenges previous reports regarding the occurrence and extent of dog admixture 
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in dingoes, as our ancestry analyses show that previous assessments severely over- 
estimate the degree of domestic dog admixture in dingo populations, particularly in 
Handling Editor: Josephine Pemberton south-eastern Australia. These findings strongly support the use of genome-wide 
SNP genotyping as a refined method for wildlife managers and policymakers to assess 


and inform dingo management policy and legislation moving forwards. 
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1 | INTRODUCTION (Allendorf et al., 2001; Todesco et al., 2016; Wayne & Shaffer, 2016). 


However, genetic admixture can also provide species with the 
Hybridisation between two populations or species leads to genetic mechanisms and genetic diversity to adapt to challenges, such as 
admixture, which alters the genetic identity of the species or popu- new or changing environments, competition or disease (Adavoudi 
& Pilot, 2022; Hamilton & Miller, 2016; Harrison & Larson, 2014). 


In many cases, genetic admixture is neutral; however, outbreeding 


lation and can result in a new hybrid species. In some cases, genetic 


dilution can lead to extinction of the species or populations involved 
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depression can result in hybrids with maladaptive genic combi- 
nations (Adavoudi & Pilot, 2022; Allendorf et al., 2001; Schmickl 
etal., 2017). 

Hybridisation and genetic admixture are of growing concern 
to conservation scientists particularly as hybridisation occurs 
more frequently between species or populations experiencing cli- 
mate change, urbanisation and/or culling (Bohling & Waits, 2015; 
Canestrelli et al., 2017; Rutledge et al., 2012). Hybridisation between 
species can lead to a hybrid swarm or unimodal hybrid zone when 
first-generation hybrids are fertile and there are no barriers to re- 
production or mating preferences (Allendorf et al., 2001; Fitzpatrick 
et al., 2015). Hybrid swarms consist of predominately intermediate 
hybrids rather than backcross progeny, or pure parental populations. 
When hybridisation occurs, but there are pre-mating barriers or as- 
sortative mating preferences, a bimodal hybrid zone may be formed, 
characterised by a population where parental forms predominate 
and intermediate hybrids are rare (Jiggins & Mallet, 2000; McFarlane 
& Pemberton, 2019). Identifying both the occurrences and reper- 
cussions of admixture can inform wildlife management and conser- 
vation policies for at-risk or keystone species (Allendorf et al., 2001; 
Fitzpatrick et al., 2015; vonHoldt et al., 2018). 

To accurately detect hybridisation and subsequent admixture, 
the genetic identity and diversity of each contributing species must 
be well described (McFarlane & Pemberton, 2019). Canis species 
(ssp.) are particularly vulnerable to hybridisation, particularly with 
domestic dogs, due to conserved genetic compatibility and repro- 
ductive behaviours, overlapping habitat and the population den- 
sity of domestic dogs (Bohling, 2016; Bohling & Waits, 2015; Lord 
et al., 2013; Macdonald et al., 2019; vonHoldt & Aardema, 2020). 
Indeed, hybridisation with free-ranging domestic dogs is a threat 
to wild canid populations in many geographic regions (Doherty 
et al., 2017; Hughes & Macdonald, 2013; Rahaman, 2017; Sykes 
et al., 2020). In this study, we sought to understand population 
structure and admixture patterns in Australian wild canids using 
genome-wide genetic markers. We hypothesise that undescribed 
population structure has biased previous admixture assessments in 
Australian wild canids, yielding significant implications for manage- 
ment and conservation (Cairns, Crowther, et al., 2021). 

The Canis (spp.) present in Australia includes dingoes, as well as 
pet and feral domestic dogs. Dingoes are an early lineage of dog, 
divergent from both modern domestic dogs and the free-ranging do- 
mestic dogs observed in Asia (Bergstr6m et al., 2020; Cairns, 2021; 
Field et al., 2022; Surbakti et al., 2020). Phylogeny and ancient DNA 
studies estimate that the dingo lineage diverged from other early 
dogs 8000-11,000years before present (YBP), prior to the rise 
of human agriculture and intensive artificial selection (Bergstr6m 
et al., 2020; Cairns & Wilton, 2016; Zhang et al., 2020). 

The management, conservation, nomenclature and taxonomy 
of dingoes are controversial, in large part due to a lack of infor- 
mation regarding a uniform definition of dingoes and the extent 
of dingoxdog hybridisation (Cairns, 2021; Crowther et al., 2014; 
Jackson et al., 2017; Smith et al., 2019). At a basic level, there are 


still key knowledge gaps concerning the genomic, morphological, 


physiological and ecological differences between dingoes and do- 
mestic dogs (Field et al., 2022; Smith et al., 2019; Eeden et al., 2018). 
In Australia, the term wild dog is often used to refer to dingoes, 
free-ranging domestic dogs and their hybrids (Crowther et al., 2014; 
Fleming et al., 2001). We hereafter use the term dingo to refer to 
pure dingoes; dingo x dog hybrids to refer to the first-generation off- 
spring of a dingo and dog; and domestic dog to refer to free-ranging 
or pet domestic dogs. We use the terms dingo backcross and dog 
backcross to refer to admixed animals with predominately dingo or 
domestic dog ancestry. Conservation policies aim to keep dingoes 
free of domestic dog admixture and minimise the impact of dingoes 
on livestock producers. While the dingo is legislatively classified 
as native wildlife, they are only protected within National Parks of 
some states and there are efforts to suppress or eradicate their pop- 
ulations across most of the Australian mainland (Bulte et al., 2003; 
Corbett, 2001; Fleming et al., 2001; Philip, 2019; Smith, 2015). 

Contemporary and historical domestic dog admixture, partic- 
ularly in south-eastern Australia, has raised significant concern 
that dingoes will become extinct through genetic dilution (Cairns, 
Crowther, et al., 2021; Claridge et al., 2014; Major, 2009; Stephens 
et al., 2015). The broad geographic range of dingoes, historical time 
period during which they and dogs have been in contact and intru- 
sive management practices that may increase dog x dingo hybridisa- 
tion are all factors thought to affect the genetic integrity of dingo 
populations (Cairns, 2021; Cairns, Crowther, et al., 2021; Stephens 
et al., 2015). 

Historically, skull morphology and pelage were used to discrimi- 
nate dingoes from dogs and identify dingo x dog hybrids (Newsome 
et al., 1980; Newsome & Corbett, 1985). However in 1999, this 
changed with the first application of DNA-based testing (Wilton 
et al., 1999), and a major continental survey (n= 3000) of dingo x dog 
hybridisation across Australia was carried out using a microsatellite 
assay comprised of 23 markers (Stephens et al., 2015). Stephens 
et al. (2015) concluded that in south-eastern Australia, <1% of the 
wild canid population were pure dingoes (Stephens et al., 2015). 
However, subsequent microsatellite DNA surveys, incorporating the 
data set of Stephens et al. (2015) and an additional 1364 samples 
from south-eastern Australia, found that most dingoes in northern, 
western and central Australia were comparatively pure (68%-98%) 
while in south-eastern Australia, only 18%-41% of animals were 
pure dingoes (Cairns et al., 2019; Cairns, Crowther, et al., 2021). 
Importantly, both Stephens et al. (2015) and Cairns, Crowther, 
et al. (2021) showed that the prevalence of free-ranging domestic 
dogs in Australia was very low with less than 1% of samples identi- 
fied as feral dogs. 

The reported low occurrence of pure dingoes in southern and 
eastern Australia has been a key factor driving government policy 
and practice (Allen et al., 2017; Bird & Bowman, 2016; DEPI, 2013; 
Major, 2009; NWDAP, 2020). Government policy and communica- 
tions widely use the term ‘wild dog’ rather than dingo to empha- 
sise the mixed ancestry of wild dingo populations and the presence 
of feral dogs (Letnic et al., 2012). Many Australians are unaware 


that dingoes are referred to as ‘wild dogs’ and as such targeted for 
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killing during pest management programmes on public lands (Eeden 
et al., 2018, 2020). The purportedly low occurrence of pure dingoes 
in south-eastern Australia led to Victoria listing dingoes as a threat- 
ened species, and New South Wales listing ‘hybridisation by feral 
domestic dogs’ as a key threatening process to dingoes. However, 
those and similar decisions were made based on limited resolu- 
tion microsatellite DNA testing. More advanced approaches using 
higher density genomic data have been used extensively to address 
similar questions in other species, but had not been applied to the 
problem of hybridisation between dingoes and dogs (McFarlane & 
Pemberton, 2019; Oliveira et al., 2015; Pilot et al., 2021; Stroupe 
et al., 2022; Szatmari et al., 2021), as we do here. 

Assessments of dingoxdog hybridisation are confounded by a 
reliance on (1) predefined and geographically restricted reference 
populations and (2) use of low numbers of variable markers. Recent 
studies utilising predominately Y chromosome or mitochondrial 
markers and small sample sizes (n<130) observed a broad pattern 
of biogeographic subdivision between dingoes in south-eastern 
Australia from those in the rest of Australia (Cairns et al., 2017, 2018; 
Cairns & Wilton, 2016; Greig et al., 2018; Zhang et al., 2020). Similar 
patterns of west to east subdivision have been observed in dingoes 
using three-dimensional cranial morphology (Koungoulos, 2020). 
However, persistent sampling gaps and extensive variation in dingo 
pelage, size and body shape suggest that additional uncharacterised 
biogeographic subdivisions exist. 

In this study, we sought to better characterise population struc- 
ture in Australian dingoes, as well as examine admixture between 
the dingo and domestic dog using genome-wide genotyping ap- 
proaches. Using over 195,474 SNPs on a commercial microarray, we 
provide compelling evidence that domestic dog admixture is uncom- 
mon in contemporary dingo populations. We further demonstrate 
that SNP-based technologies are a comparatively cheaper and more 
rigorous method for determining the composition of dingo popula- 
tions across Australia. In doing so, we provide the first detailed study 
of geographic population structure in dingoes across all regions of 


Australia. 


2 | MATERIALS AND METHODS 
2.1 | Canid sampling 


We collected tissue, blood or buccal samples from 307 wild and 84 
captive dingoes from locations across Australia (Figure 1, Dryad 
dataset: [https://doi.org/10.5061/dryad.vq83bk3ww)]). A majority 
of the wild dingoes were sampled as part of programmes to suppress 
dingo populations. Additional samples were provided by private 
landholders, government officers, found as roadkill or during car- 
cass disposal or sampled from puppies found in the wild. Unrelated 
captive dingo samples were contributed by dingo sanctuaries, zoos 
and private citizens. We included captive dingo samples to serve 
as controls in a comparison to wild dingoes, due to their known 


parentage. A subset of the wild and captive samples were from 


existing scientific collections and had been previously tested using 
a 23-marker microsatellite-based test developed for distinguishing 
dingoes from dingoxdog hybrids (Cairns, Crowther, et al., 2021; 
Stephens et al., 2015). We also collected blood or buccal samples 
from 36 Australian domestic pet dogs and 111 North American 
domestic pet dogs. We chose dog samples from working, herd- 
ing or mixed breeds as they represent the populations most likely 
to have hybridised with Australian dingoes. A set of five unrelated 
Vietnamese putative primitive breed samples were sourced from a 
private owner in H6 Chi Minh City, Vietnam. DNA was extracted 
from blood, tissue or buccal samples using either (1) Qiagen DNeasy 
Blood and Tissue (Qiagen Sciences) or (2) Monarch Genomic DNA 
Purification kits (New England Biolabs) or (3) a modified proteinase 


K, phenol-chloroform extraction method (Bell et al., 1981). 


2.2 | Axiom canine genotyping 


After DNA quantification, the dingo, Australian dog and Vietnamese 
dog samples were genotyped on the Axiom Canine HD Genotyping 
array (Thermo Fisher Scientific Inc.) at the Ramaciotti Centre 
for Genomics (University of New South Wales). The 111 North 
American domestic pet dogs were genotyped on Axiom Canine 
Genotyping Array Set A and Set B (Thermo Fisher Scientific Inc.) at 
the Thermo Scientific Microarray Research Services Laboratory. We 
chose to use the >700,000 marker Axiom Canine HD Genotyping 
array (Thermo Fisher Scientific Inc.) rather than the >170,000 
marker Illumina CanineHD Whole-Genome Genotyping BeadChip 
(IIlumina Inc.) due to the increased number of informative markers 
(>130,000) contained in the Axiom Canine arrays. Genotypes were 
called using Axiom Analysis Suite (Thermo Fisher Scientific Inc.) and 
quality control (QC) filtering was conducted using Plink v1.9 (Purcell 
et al., 2007). We trimmed the North American domestic dog data to 
include only the SNPs from the Axiom Canine HD Genotyping array 
(Thermo Fisher Scientific Inc.), as the Axiom Canine Genotyping Set 
A and Set B arrays (Thermo Fisher Scientific Inc.) contain probes for 
a larger set of SNP markers and were used in the development of the 
Axiom Canine HD Genotyping array (Thermo Fisher Scientific Inc.). 
We excluded individuals or sites with 210% missing data using the 
commands --geno 0.1 and -- missing. We did not perform minor allele 
frequency filtering because it has a negligible impact on ascertain- 
ment bias and may negatively impact genetic diversity or population 
structure assessments. Instead, we employed linkage disequilibrium 
(LD) filtering (Malomane et al., 2018), by removing regions in high LD 
using the command --indep-pairwise 50 5 0.6 in Plink v1.9. 


2.3 | Dingo x dog ancestry and admixture analysis 


We performed genetic ancestry analyses on the complete data set of 
307 wild dingoes, 84 captive dingoes and 152 domestic breed dogs 
using FastStructure (Raj et al., 2014). We employed FastStructure, 
a Bayesian population clustering method similar to STRUCTURE 
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FIGURE 1 Map drawn in QGIS depicting the sampling distribution of 307 wild dingoes across Australia. Dashed lines depict the 


geographical sampling regions defined as Western Australia (1), Central Australia (2), Northern Australia (3), Eastern Australia (4), Southern 
Australia (5) and the Big Desert (6). The Australian States are Western Australia (WA), Southern Australia (SA), the Northern Territory (NT), 
Queensland (QLD), New South Wales (NSW), Victoria (VIC), Tasmania (TAS) and the Australian Capital Territory (ACT). The solid black lines 
indicate the position of Government maintained dingo exclusion fences. 


(Pritchard et al., 2000), due to its capacity to handle large genomic 
data sets and integrated method of choosing model complexity (K- 
values). Wild dingoes belonged to one of six populations according to 
geographic sampling location: Western Australia, Central Australia, 
Northern Australia, Eastern Australia, Southern Australia and the Big 
Desert (Figure 1). Captive-bred dingoes were identified as such for the 
analysis and domestic dogs were identified as either a specific breed or 
a mixed breed group based on owner reporting. The specific domestic 
dog breed groups were Vietnamese dog (VIET), Mixed breed (MIXED), 
Kelpie (KELPIE), Australian Cattle Dog (ACD), Australian Stumpy Tail 
Cattle Dog (ASTCD), Border Collie (BORD), Golden Retriever (GOLD), 
German Shepherd (GSD), Labrador (LAB), Patagonian sheepdog 
(PGOD), Tervuren Shepherd (TURV) and Belgian Shepherd (BELS). 
FastStructure analyses were ‘unsupervised’ meaning that we did not 
use a-priori defined populations to inform the Bayesian modelling and 
admixture detection (Raj et al., 2014). 


FastStructure analyses were run with the StructureThreader pipe- 
line (Pina-Martins et al., 2017), modelling K=1-20 with a random seed 
and simple prior. A simple prior was suitable for our study because we 
are primarily interested in strong population structure between din- 
goes and domestic dogs rather than substructure between domestic 
dog breeds. FastStructure reports the optimal K values to maximise 
marginal likelihood and explain structure in the data (Raj et al., 2014). 
Ten replicates were completed for each modelling scenario that best 
fit the data (K=8-10), calculating the average proportion of dingo an- 
cestry for each individual along with the standard error. Confidence 
intervals were calculated by multiplying standard error values by 1.96 
(McFarlane et al. 2020). The proportion of estimated dingo ancestry 
for each individual with 95% confidence intervals was plotted in R 
using the ggplot package (Wickham, 2016). An ancestry component 
bar chart based on the average FastStructure q-values was created in 
R using the BITE package (Milanesi et al., 2017). 
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FastStructure ancestry component q-values were used to cate- 
gorise each animal as a dingo, historical or recent dingo backcross, 
dingo x dog hybrid, recent or historical dog backcross or domestic dog 
(Table 1). Rather than using a strict threshold, we assigned animals 
to a specific category using the calculated 95% confidence intervals 
(McFarlane & Pemberton, 2019). We considered that the occurrence 
of different backcross classes in wild dingo populations may inform 
conservation and management policy directions, with some previous 
studies considering animals with >93% dingo ancestry as pure (Allen 
et al., 2017). Therefore, we categorised backcrosses as being either: 
(1) historical backcrosses (>93% parental ancestry) that were likely 
to be third generation or higher or (2) recent backcrosses (55%-93% 
parental ancestry) which described animals that were likely to have 
had a dingo x dog hybrid ancestor in the last one to three generations 
(Dziech, 2021; Harmoinen et al., 2021; McFarlane & Pemberton, 2019). 
A map depicting the FastStructure ancestry assignment of each wild 
origin sample was created using the QGIS v3.12 (QGIS, 2020). 

Principal coordinates analyses (PCoA) were carried out using the R 
package SambakR function ‘find_structure’. In the SambaR analyses, pop- 
ulations were defined according to the FastStructure results with ances- 
try categorised as in Table 1. However, Vietnamese dogs were allocated 
to their own category, given that dogs from South Asia and dingoes 
may share ancestry (Cairns, 2021). F., values for animals categorised as 
dingoes, (historical and recent) dingo backcrosses, dingo x dog hybrids, 
(historical and recent) dog backcrosses, domestic dogs and Vietnamese 


dogs were calculated in SambaR using the ‘calcdistance’ function. 


2.4 | Dingo population structure and biogeography 


To explore patterns of biogeographic variation in dingoes across 
Australia, the SambaR ‘find_structure’ function was used to run 
PCoA using only animals with an ancestry category of dingo or his- 
torical dingo backcross, excluding all animals with greater than 7% 
ancestry deriving from domestic dogs. Analyses were run with and 


without dingoes of captive origin. A map of wild origin dingo samples 
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and their primary FastStructure population cluster ancestry identity 
was plotted using QGIS v3.12 (QGIS, 2020). 

Inbreeding coefficients (F,.,) were calculated in Plink v 1.9 using 
the command --het based on the equation F,,,.=(O - E)/(N - E) where 
O is the number of observed homozygous markers of the individ- 
ual, E is the expected number of homozygous markers under the 
Hardy-Weinberg equilibrium calculated from the allele frequencies 
estimated on the dingo data set and N is the total number of markers 
(Gazal et al., 2014; Purcell et al., 2007). F,., is an individual-based 
inbreeding estimator allowing us to compare differences in homo- 
zygous excess resulting from inbreeding between the identified 
FastStructure dingo populations, historical dingo backcrosses, re- 
cent dingo backcrosses and dingoxdog hybrids. All samples iden- 
tified as domestic dogs, Vietnamese dogs or (recent or historical) 
dog backcrosses were excluded from this analysis. Animals were as- 
signed to subpopulations based on the primary FastStructure dingo 
population cluster and boxplots visualising the range of computed 
F. 4 Values in the different dingo subpopulations were created in R 


using the ggplot package (Wickham, 2016). 


2.5 | Comparisons between microsatellite and SNP 
ancestry estimates 


We compared the modelled ancestry q-values for a set of 113 ani- 
mals that were tested using both the 195,434-marker Axiom SNP 
array method and the 23-marker microsatellite method of Cairns, 
Crowther, et al. (2021) and Stephens et al. (2015). We calculated the 
difference between our modelled SNP ancestry q-values and previ- 
ously reported values based on microsatellite DNA testing (Cairns, 
Crowther, et al., 2021; Stephens et al., 2015). We then calculated 
the average difference between SNP and microsatellite q-values for 
individuals within the five dingo subpopulations. We also assessed 
the number of individuals which were categorised by microsatellite 
DNA testing as recent backcrosses or dingo x dog hybrids, but were 


pure or historical backcrosses based on SNP analysis. 


TABLE 1 Genetic thresholds used for the categorisation of dingo ancestry during microsatellite* or SNP (this paper) DNA testing 


methods. 


Microsatellite q-value 


Microsatellite category* threshold? 
Dingo 1 (Dingo with no dog ancestry) >.90 
Dingo 2 (likely dingo with no dog ancestry) .80-.90 
Dingo with dog ancestry 1 (>75% dingo) .70-.80 
Dingo with dog ancestry 2 (65%-75% dingo) .60-.70 
Dingo with dog ancestry 3 (55%-64% dingo) .55-.60 
Dingo x dog hybrid (likely F1 hybrid) 45-.55 
Domestic dog with dingo ancestry (<50% dingo) 45-.25 
Domestic dog .25-.0 


SNP dingo q-value 


Genomic (SNP) category threshold 
Dingo >99 
Historical dingo backcross 93-99 
Recent dingo backcross .55-.93 
Dingo x dog hybrid .45-.55 
Recent dog backcross .07-.45 
Historical dog backcross .07-.01 
Domestic Dog <.01 


°As per Stephens et al. (2015) and Cairns, Crowther, et al. (2021), Cairns, Newman, et al. (2021). 
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3 | RESULTS 
3.1 | Axiom canine genotyping 


Genotyping of the 543 canid samples on the Axiom Canine 
Genotyping arrays (Thermo Fisher Scientific Inc.) resulted in raw 
genotypes at 627,324 variable loci (Cairns et al., 2023). After QC and 
LD filtering using Plink v1.9, a data set of 195,474 SNPs was available 
for dingox dog admixture and biogeographic clustering analysis in 
FastStructure and SambaR (Cairns et al., 2023). 


3.2 | Dingo x dog ancestry and admixture analysis 


FastStructure ancestry modelling identified K=8 as the best fit for 
the data while K=10 maximised model complexity. We observed 
the presence of five biogeographic dingo population clusters which 
we termed: Big Desert, West, East, South and Captive (Figure 2). 
Dingoes from central and northern Australia carried a mixture of 
West and East population ancestry and did not form distinct popula- 
tion clusters. Some dingoes from eastern Australia carried a mixture 
of East and West population ancestry. Domestic dog breeds did not 
cluster discretely from one another; however, three main clusters 
were observed: (1) Border Collie, (2) Belgian Shepherd/Turvian 
Shepherd and (3) all other dogs. The Vietnamese primitive breed 
dogs did not form their own population cluster but showed ancestry 
from both domestic dogs and dingoes. 

Each sample was assigned using the Table 1 SNP testing criteria 
to an ancestry category after assessing calculated 95% confidence 
interval dingo q-values (Figure 3). In circumstances where the 95% 
confidence interval spanned a category threshold, we used the 
lower bound of the 95% confidence interval to classify individuals. 
This allowed us to conservatively assign animals and minimise the 
number of Type II errors where a backcross animal is assigned as 
a parental type. Within our data set of 307 putative dingoes from 
across Australia, 70.0% were pure dingoes; 15.0% were histori- 
cal dingo backcrosses; and 15.0% were recent dingo backcrosses 
(Table 2). No wild dingo x dog hybrids or feral domestic dogs were 
identified. All seven of the dingo x dog hybrids in our data set were 
born and kept in captivity. Recent dingo backcrosses were more 
common in QLD (26.5%) and NSW (20.4%), while there were few 
recent dingo backcrosses observed in VIC (6.5%). Indeed, the four 
recent dingo backcrosses from VIC all carried >85% dingo ancestry. 
There was minimal evidence of dog ancestry in the dingoes sampled 
from SA or WA with only a single recent dingo backcross observed. 
From the NT, there were two recent dingo backcrosses collected 
near towns, and the rest (n=14) were pure dingoes or historical 
dingo backcrosses. 

PcoA provided additional evidence of differentiation between 
dingoes and domestic dogs, while highlighting the limited amount 
of dog ancestry in dingoes (Figure 4). Admixed dingoes clustered 
closer to pure dingoes than to domestic dogs. Dingo x dog hybrids 


were distributed roughly equidistant between the pure dingoes and 


domestic dog clusters. Interestingly, the Vietnamese dogs behaved 
similar to dingo x dog hybrids, corroborating the FastStructure anal- 
ysis (Figures 2 and 3), revealing both dingo and dog ancestry. K=10 
modelling suggests there may be further regional substructure 
within dingo populations, particularly the East population cluster, as 
we observed five animals from the NSW Southern Highlands near 
Sydney forming their own cluster. 

F,, values indicate that dingoes are distinct from domestic 
dogs, F,,=0.288; domestic dog backcrosses, F,,=0.293-0.324; 
and Vietnamese dogs, F,,=0.301 (Table 3). There was a moderate 
level of differentiation between dingoes and dingoxdog hybrids 
with F,,=0.126. Dingo backcrosses while not distinct from pure 
dingoes (F,,=0.008-0.028) were distinct from both domestic dogs 
(F,;=0.157-0.259) and Vietnamese dogs (F,,>0.191). Dingo back- 
crosses were only mildly differentiated (F,,=0.039-0.088) from 
dingoxdog hybrids (Table 3). Calculated inbreeding coefficient 
(Fret) values indicate that dingoes, both pure and backcrossed, had 
elevated homozygosity levels compared to dingoxdog hybrids 
(Figure 5). Captive dingoes had similar F,., values to wild dingo 


populations. 


3.3. | Dingo population structure and biogeography 


Mapping of the FastStructure modelled primary population cluster 
identity for the 307 wild and 84 captive dingoes depicts a pattern of 
subdivision between the five population clusters: West, South, East, 
Big Desert and Captive (Figure 6). Dingoes in northern and central 
Australia carry a mixture of East and West population ancestry with 
the proportion of West dingo ancestry increasing north along the 
QLD coast or inland away from the east coast. Dingoes that originate 
from captive breeding largely form their own population, but there 
were three instances of putatively wild dingoes carrying the captive 
population cluster ancestry. 

The population substructure observed in FastStructure is 
corroborated by the PcoA with four wild dingo populations: Big 
Desert, South, East and West (Figure 7). The Captive dingo popu- 
lation forms its own cluster. The large variation within population 
clusters (Figure 7) may corroborate the observation of additional 
substructure in the K=10 FastStructure modelling (Figure 2). 

Calculated F,, values indicate that there has been some gene 
flow between the West, East and South dingo populations (Table 4). 
Indeed, dingoes in northern and central Australia carry mixed East 
and West ancestry while some South dingoes carry some West 
dingo ancestry (Figure 2). The Big Desert population is the most dis- 
tinct dingo population with limited evidence of gene flow between 
it and any of the other dingo populations. The Captive population 
cluster also had elevated F,, values suggesting limited geneflow be- 
tween captive and wild dingo populations. 

Comparison of median F,., values between the dingo population 
clusters indicates that all dingoes have higher homozygosity levels 
than recent dingo backcrosses and dingoxdog hybrids (Figure 5), 


which is consistent with whole genome sequencing data that also 
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FIGURE 2 FastStructure modelling of 543 canids and 195,474 genomic SNPs for K=8-10. Samples grouped according to geographic 
region or owner reported breed. The dingo populations are defined according to geographical sampling region: West Australia (WEST), 
Central Australia (CENTRAL), Northern Australia (NORTH), Eastern Australia (EAST), Southern Australia (SOUTH) and Big Desert 
(BIGDESERT) as per Figure 1. Dingoes from captive-bred origin were grouped together as CAPTIVE. Domestic dog samples were defined 
as per owner-described breed groups including Vietnamese dog (VIET), Mixed breed (MIXED), Kelpie (KELPIE), Australian Cattle Dog (ACD), 
Australian Stumpy Tail Cattle Dog (ASTCD), Border Collie (BORD), Golden Retriever (GOLD), German Shepherd (GSD), Labrador Retriever 
(LAB), Patagonian sheepdog (PGOD), Tervuren Shepherd (TURV) and Belgian Shepherd (BELS). 


reported higher levels of homozygosity in dingoes compared to 
domestic dogs (Kumar et al., 2023; Zhang et al., 2020). The Big 
Desert population has the highest median F,., values, nearly 30% 


higher than the other dingo populations, indicating the population 


is extremely homozygous and likely inbred. We did not compare 
homozygosity between dingoes and domestic dogs due to potential 
ascertainment bias associated with the design of the Axiom Canine 


HD Genotyping array (Thermo Fisher Scientific Inc.). 
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FIGURE 3 Average Dingo ancestry 
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TABLE 2 Observation of dingo ancestry categories across Australia in 307 wild and 84 captive dingoes based on SNP genomic testing. 


Wild Captive 
WA QLD NSW? vIC NT SA Australia Australia 
% (n) % (n) % (n) % (n) % (n) % (n) % (n) % (n) 
Dingo (>99% dingo) 83.3 (37) 54.4 (37) 62.1 (64) 87.1 (54) 68.8 (11) 87.5 (14) 70.0 (215) 63.1 (53) 
Historical dingo backcross 14.3 (6) 19.1 (13) 17.5 (18) 6.5 (4) 18.8 (3) AS) (P) 15.0 (46) 7.1 (6) 
(>93% dingo) 
Recent dingo backcross 2.4 (1) 26.5 (18) 20.4 (21) 6.5 (4) 12.5 (2) - 15.0 (46) 13.1 (11) 
(55%-93% dingo) 
Dingo-dog hybrid (45%- - - - - - - - 8.3 (7) 
55% dingo) 
Recent dog backcross a 7 - - - = 7 8.3 (7) 
(<55%-93% dog) 
Historical dog backcross - - - - - - - - 
(>93% dog) 
Domestic dog (>99% dog) - - - - - - - - 
Total (n) 42 68 103 62 16 16 307 84 


"Includes one sample from the Australian Capital Territory (ACT). 


3.4 | Comparisons between microsatellite and SNP 
ancestry estimates 


Nuclear SNP-based approaches are now the preferred method for an- 
cestry assessment of mammals throughout the tree of life (McFarlane 
et al., 2020; Steyer et al., 2018; Stroupe et al., 2022; Szatmari 
et al., 2021; Vali et al., 2010; Zimmerman et al., 2020). However, 
as microsatellites have historically been the method of choice, we 
conducted a direct comparison between SNP and _ microsatellite- 


calculated dingo ancestry, observing that the proportion of dingo 


ancestry and occurrence of dingo x dog hybridisation across Australia 
differs depending on the method used (Figure 8, Tables 2 and 5). The 
microsatellite DNA testing method severely overestimates the pres- 
ence of dog ancestry in dingoes compared to the SNP DNA testing 
method (Figure 9). For animals from the Captive, Big Desert, East or 
South population clusters, microsatellite analysis results in the as- 
signment of animals to the incorrect category 28%-44% of the time 
(Table 5). While microsatellite DNA testing demonstrates widespread 
dog admixture in south-eastern Australia, SNP data show limited do- 


mestic dog ancestry in the wild dingo population. 


asuaor’] SUOWIUIOD aANearD a[qvorfdde oy) Aq pousoAos are sajotue YO ‘asn Jo sapns Joy Aresqr’y] ouTpUC A2|I AA UO (SUONIpUOD-puL-suLIa}/WIOD KoTTM Aresgyauruo//:sdyyy) SUOTIPUOD puL SULIAL, aI 99g ‘[PZOZ/10/9T] UO Aresqry auTTUG AIIM ‘SGOT OOW/T TT LOL opstutos: Kopi: Arerqryaurfuoy/:sdyyy wos papeopumog ‘ST ‘EZOT ‘XP6TS9ET 


CAIRNS eT AL. 4141 
MOLECULAR ECOLOGY gaAY@| Ley 
PC1 vs PC2 PCl1 vs PC3 PC2 vs PC3 
@ Dingo . e ad 
o| M Historical Dingo backcross ? Pal 3 val ¢ 3 
= ~~ Recent Dingo backcross S oot o| 
©]  m Dingo x dog hybrid @ 
™ Vietnamese dog e e e e@ ae 
@ Recent dog backcross e e e 
«| @ Historical dog back _~ e ~ 
a ° (S o “eS 
Nae e S = ai = Co e 
N e* e — @ VY e 
[e) 8 t | e isa) foal 
el peer f° = e r . 
o 
” T T 8 
Co 
: é 4, 
e e 
0.10 — -0.05 0.00 0.05 0.10 0.15 0.10 0.05 0.00 0.05 0.10 0.15 0.05 0.00 0.05 0.10 
PCI (41.9%) PC1 (41.9%) PC2 (5%) 
FIGURE 4 Principal coordinates analysis (PcoA) based on Nei distances for 543 canids based on 195,474 genomic SNPs constructed in 
SambaR. Samples are coloured based on FastStructure ancestry classification. 
TABLE 3 Pairwise F,, values calculated for the five dingo versus dog ancestry categories* 
Historical dingo Recent dingo Dingo x dog Recent dog Historical dog Domestic Vietnamese 
backcross backcross hybrid backcross backcross dog dog 
Dingo .008 .028 126 .293 324 .288 301 
Historical dingo - .012 .088 LEO SD) a2 P25 
backcross 
Recent dingo - - .039 {157 181 171 191 
backcross 
Dingo x dog hybrid = = = 046 .053 065 AS} 
Recent dog - - - - .012 .029 116 
backcross 
Historical dog - - - - - .012 oll 
backcross 
Domestic dog - - - - - - 143 


*Fot values calculated in SambaR using data set of 195,474 SNPs. 
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FIGURE 5 F,., (inbreeding coefficient) a measure of homozygosity calculated in Plink v1.9 from 195,474 SNPs, as observed in the five 
distinct dingo population clusters, Big Desert, West, East, South and Captive. F,,, was also calculated for historical dingo backcrosses, recent 
dingo backcrosses and dingo ~ dog hybrids. Vietnamese dogs, domestic dogs, historical dog backcrosses and recent dog backcrosses were 
excluded from this analysis due to possible ascertainment bias. 
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4 | DISCUSSION 


Ancestry modelling, conducted using high-density genome-wide 
SNPs, reveals two key findings: (1) there are at least four wild dingo 
populations and a separate captive dingo population in Australia; (2) 
the presence of dog ancestry in wild dingoes is much less common 
than previously hypothesised by microsatellite DNA testing (Cairns, 
Crowther, et al., 2021; Stephens et al., 2015). These data provide 
much needed clarity regarding the identity of dingoes across Australia, 
which can inform policy debate regarding how these animals are man- 
aged. Furthermore, these findings indicate that genetic surveys of wild 
dingo populations using genome-wide technologies are important for 
understanding dingoxdog ancestry across Australia. Future use of 
genome-wide technology, particularly restriction-site associated DNA 
sequencing (RADseq) or whole genome sequencing, will permit simul- 
taneous analysis of domestic dog admixture, demography, homozy- 
gosity and genomic diversity in dingo populations while also enabling 
comparisons to domestic dogs or wolves. Finally, these data argue for 


significant changes in current genetic testing methods. 


ecoeoe t—t—“‘“‘ ‘“ ‘ROO 
Captive 


4.1 | Biogeography and regional variation 
in the dingo 


We observed the presence of four distinct geographically distributed 
dingo populations: West, East, South, Big Desert (Figures 2, 6 and 7). 
Captive dingoes formed their own population, distinct from wild dingo 
populations (Figures 1 and 6). The Big Desert dingo population was the 
most strongly differentiated, with F,, values twice that of other popu- 
lations (Table 4). While there are several phenotypically described va- 
rieties of dingo, it is unclear if these align with the genetic populations 
(Corbett, 2001; Walters, 1995). Our biogeographic findings corrobo- 
rate and extend the observation of geographic subdivision in dingoes 
from skull morphometrics, mitochondrial DNA, Y chromosome and 
nuclear DNA data (Cairns et al., 2017, 2018; Cairns & Wilton, 2016; 
Greig et al., 2018; Koungoulos, 2020; Zhang et al., 2020). We hypoth- 
esise that there is additional regional substructure within the four 
geographic dingo populations, with preliminary evidence of a region- 
ally distinct subtype of the East dingo population around the NSW 
Southern Highlands, southwest of Sydney (Figure 2). 


Majority dingo ecotype identity 
Big Desert © 


West e 
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South ce) 
Captive e 
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FIGURE 6 Map drawn in QGIS of the primary dingo population identity of 320 dingo and historical dingo backcross samples from the 
wild or captivity based on FastStructure modelling of 195,474 SNPs. Samples are coloured according to primary dingo population cluster 
identity as assigned by FastStructure. The solid black lines indicate the position of Government maintained dingo exclusion fences. 
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FIGURE 7 Principal coordinates analysis (PcoA) of Nei distances carried out in SambaR of 320 dingo and historical dingo backcross 
samples with less than 7% dog ancestry based on 195,474 genomic SNPs. Samples are coloured according to primary dingo population 


cluster identity as assigned by FastStructure. 


TABLE 4 Pairwise F,, values between the five dingo 
populations® 


West East South Captive 
Big Desert .276 .267 .29 343 
West = .056 .084 .098 
East - - .063 096 
South = = = odl 


*Fot values calculated in SambaR using data set of 195,474 SNPs. 


Geneflow between the distinct wild dingo populations may re- 
flect differing dispersal and demographic patterns across Australia. 
Our analyses observed geneflow between the West and East dingo 
populations in central and northern Australia and in southern 
Australia between the South and East dingo populations (Figure 2, 
Table 4). However, there was limited evidence of geneflow into or out 
of the Big Desert dingo population (Figure 2). The pattern of strong 
genetic differentiation (F,,=0.267-0.343) and limited geneflow 
between the Big Desert and other dingo populations likely reflects 
the influence of historical and ongoing lethal control programmes in 
the surrounding area. Lack of knowledge concerning historical and 
contemporary biogeographic patterns in dingoes warrants further 
investigation using DNA samples from ancient or colonial-era indi- 
viduals, and comparison to Oceanic canids such as the New Guinea 
Singing Dog (Surbakti et al., 2020). 

Dingoes from all four of the wild populations were observed 
to have increased levels of homozygosity (median F,,>0.3) rela- 
tive to recent dingo backcrosses and dingo x domestic dog hybrids 
(Figure 5). We observed elevated homozygosity levels (median 
F.4.> 0.7) in the Big Desert compared to other dingo populations, 
suggesting the population is experiencing inbreeding (Figure 5). A 
recent microsatellite-based study found the Big Desert dingo pop- 
ulation to have genetic diversity levels up to 50% lower than other 
dingo populations (Stephens et al., 2022). While the absolute value 


of homozygosity measures calculated for dingoes in this study could 
be inflated due to ascertainment bias, whole genome data from a 
small number of captive and wild dingoes also reveal elevated homo- 
zygosity levels relative to domestic dogs (Field et al., 2022; Kumar 
et al., 2023; Zhang et al., 2020). Further research is needed to as- 
sess whether lethal management programmes have had an impact 
on homozygosity and genetic diversity levels in regional dingo pop- 
ulations, such as the Big Desert population (Hohenlohe et al., 2021; 
Quevedo et al., 2019; Weeks et al., 2016). The increased levels of 
homozygosity (median F,,,,>0.3) found in dingoes may also reflect 
ancestral bottlenecks, which could be investigated in the future 


with runs of homozygosity (ROH) analysis of genomic data. 


4.2 | Domestic dog admixture in Australia 


We herein provide the most accurate data regarding the prevalence 
of domestic dog ancestry in the dingo population, particularly in 
south-eastern Australia. Our survey of 307 wild dingoes from across 
Australia revealed that the occurrence and extent of domestic dog 
ancestry within the dingo population was limited and no feral do- 
mestic dogs were observed (Figures 2 and 3). In Victoria, where pre- 
vious studies using microsatellites suggest that the dingo population 
is small and highly admixed (Cairns, Crowther, et al., 2021; Stephens 
et al., 2015), we observed that 87.1% (n=54) of animals tested were 
pure dingoes, and 6.5% (n=4) were historical dingo backcrosses, with 
>93% dingo ancestry (Table 2). Similarly, in New South Wales and 
Queensland (n=171), where dingo x domestic dog hybridisation had 
been considered pervasive (Cairns, Crowther, et al., 2021; Stephens 
et al., 2015), we only observed two wild canids with <70% dingo 
ancestry. We found that most of the sampled New South Wales and 
Queensland animals were pure dingoes (n= 101) or historical dingo 
backcrosses (n=31). In the Northern Territory, South Australia and 
Western Australia, we found little evidence of domestic dog admix- 


ture in the dingo population. 
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FIGURE 8 Comparison between estimated proportion of domestic dog ancestry observed in a set of 113 wild and captive dingoes 
using both SNP and microsatellite ancestry testing methods. Samples are grouped according to primary dingo population cluster identity as 


assigned by FastStructure: Big Desert, West, East, South and Captive. 


TABLE 5 Comparison of individual SNP and microsatellite DNA testing results for a set of 113 wild and captive dingoes. 


Australia 
Number of animals tested by both methods 113 
Mean difference in SNP and microsatellite DNA a5 


modelled dingo ancestry values (q-value) 


Percent of animals where SNP and microsatellite DNA 20% 
testing disagrees on ancestry categorisation 


Our genome-wide SNP data demonstrate that dingoes and do- 
mestic dogs are genetically distinct (Figures 3 and 4). F,, values 
indicate contemporary and ongoing geneflow between admixed 
and non-admixed dingoes. Ongoing hybridisation between dingoes 
and dogs of different breeds would disrupt patterns of geographic 
population structure and decrease homozygosity. The geographic 
structure and elevated homozygosity levels observed in dingoes, 
particularly in southern and eastern Australia, provide additional 
evidence that there is limited admixture between dingoes and do- 
mestic dogs (Figures 4-6). Dingoes from southern Australia exhibit 
greater phenotypic variation than those elsewhere, particularly in 
terms of pelage (Cairns, Newman, et al., 2021; Jones, 1990, 2009; 
Newsome & Corbett, 1985). While variable pelage was previously 
hypothesised to be the result of domestic dog admixture (Newsome 
& Corbett, 1985), our finding of rare domestic dog admixture makes 
this unlikely, and suggests that diverse coat colours may repre- 
sent standing ancestral variation or perhaps local adaption (Cairns, 
Newman, et al., 2021). 

It was hypothesised that dingoxdomestic dog admixture pat- 
terns in Australia fit a unimodal hybrid swarm with extensive ge- 
netic swamping (Allen et al., 2017; Claridge et al., 2014; Daniels & 
Corbett, 2003; Stephens et al., 2015). However, our data suggest that 
the pattern of genetic admixture observed in dingoes fits a bimodal 


hybrid zone, with a majority of the population being pure dingoes 


Big desert Captive East South West 
16 3 29 13 52 
AA .18 Ale 16 09 
44% 33% 28% 31% 6% 


or dingo backcrosses with low levels of intermediate dingox dog 
hybrids (Figures 3 and 4, Table 2). The low occurrence of both first- 
generation dingo x dog hybrids and low dingo ancestry backcrosses 
further supports this hypothesis as it suggests the presence of 
behavioural barriers or assortative mating preferences (Allendorf 
et al., 2001; Fitzpatrick et al., 2015; McFarlane & Pemberton, 2019). 
It has also been suggested that lethal management strategies such 
as aerial 1080 baiting, particularly during the dingo breeding sea- 
son, increase the risk of dingo x dog hybridisation by fracturing so- 
cial structures and reducing the availability of dingo mates (Cairns 
et al., 2019; Cairns, Crowther, et al., 2021; Wallach et al., 2009). 
This might explain why dog ancestry is more commonly observed 
in dingo populations where intensive or aerial baiting strategies are 
widespread, such as NSW and QLD (Table 2). Further research is 
needed to identify whether behavioural barriers to hybridisation or 
assortative mating preferences occur in dingo populations and to 
understand the impact of lethal control on patterns of dog admix- 


ture in dingo populations. 


4.3 | Dog ancestry detection methodology 


Our results demonstrate that microsatellite-based DNA methods 


consistently estimate higher levels of dog ancestry in dingoes and 


asuaor’] SUOWIUIOD aANearD a[qvorfdde oy) Aq pousoAos are sajotue YO ‘asn Jo sapns Joy Aresqry] aulpuC Ad|IA\ UO (SUONIpUOD-puL-suLIA\/WIOD KoTTM Aresgyauruo//:sdyy) SUOTIPUOD pure SUIAL By 99g ‘[PZOZ/10/9T] UO Aresqry ouTTUG AIIM ‘8G69T'OOW/T TT LOL opsuros: Kopi: Arerqryaurfuoy/:sdyyy wos papeopumog ‘CT ‘EZOT ‘XP6TSIET 


CAIRNS eT AL. 4145 
MOLECULAR ECOLOGY Sai sate ea 
SNP Ancestry categories Microsatellite Ancestry categories 
(a) @ Dingo (b) © Dingo 1 & 2 
- © Historical Dingo backcross ° e Dingo with dog ancestry 1 
Darwin ® Recent Dingo backcross 2 ar ~ Dingo with dog ancestry 2 
| q e pee La esi * a q Dingo with dog ancestry 3 
e ® .___@ Recent Dog backcross “mae of © Dingo x Dog hybrid 
Sr. mecais Historical Dog backcross $: 3 38 © Dog 
ss Lr ©, © Domestic Dog 
at e e e 
eid | ele . 
Alice Springs @ 
e e 
e os e Soprisbane 
e 3 LA A o's 
e 
oe %e ou o 2 
Pert! a 
Adelaide © Sydney a © canberas3° Sine 
-) Canberra ‘ew ‘ y 
nel Melbourne: 
0 250 500k A a ag A 
= _—— 
(c) “Hobart (d) 0 250-500» 750-—1000km Poort 
100;— —_ = = - - 100 a | 
90 | oa 90 
80 : 80 
70 i 70 
60 60 
40 40 
30 30 
20 20 
10 10 
0 oSaeay! = GaN ee NRE eeGeey, haere oo ese Gla esrcagee 0 
WA QLD NSW&_ VIC NT SA Australia WA SA NT QLD NSW VIC ACT 
ACT 


FIGURE 9 Comparison of the occurrence and prevalence of dingo versus dog ancestry across Australia based on a 307 sample SNP 
data set with 195,474 genomic markers (A, C) and a 5036 sample data set with up to 23 microsatellite genetic markers (B, D) from Cairns, 


Crowther, et al. (2021). 


miscategorised pure dingoes as dingo backcrosses compared to 
analysis based on 195,474 SNPs (Figure 9). Interestingly, dingoes 
from the Big Desert, South, East or Captive population clusters 
were more often miscategorised by microsatellite-based methods 
than those from the West (Table 5). Our findings are consistent 
with reports from a wide range of other species demonstrating that 
genome-wide SNP analysis outperforms microsatellites in terms of 
accurate identification of ancestry and admixture (Dziech, 2021; 
Mattucci et al., 2019; McFarlane et al., 2020; Steyer et al., 2018; 
Stroupe et al., 2022; Szatmari et al., 2021; Vaha & Primmer, 2006; 
Vali et al., 2010; Zimmerman et al., 2020). 

While microsatellite-based markers remain useful for short 
evolutionary timeframes, or in taxa that are slow evolving or clonal 
(Putman & Carbone, 2014), they are much less useful for studies 
of population genetics and demography, hybrid detection or anal- 
ysis of large pedigrees (McFarlane & Pemberton, 2019; Putman & 
Carbone, 2014; Szatmari et al., 2021; Telfer et al., 2015; Tokarska 
et al., 2009; Vali et al., 2008; Thaden et al., 2017; Zimmerman 
et al., 2020). We suggest that microsatellite DNA testing under- 
estimates dingo ancestry because (1) it relies on geographically 
restricted reference populations with strict thresholds, (2) uses 
a relatively small number of genetic markers that are not strictly 
diagnostic and (3) the Bayesian modelling is typically constrained 
to scenarios with K=2. The microsatellite markers used for an- 


cestry assessment in dingoes were chosen by comparing allele 


frequencies between a small core set of mixed breed domestic 
dogs and dingoes from captivity or western and central Australia 
(Cairns, Crowther, et al., 2021; Elledge et al., 2008; Stephens 
et al., 2015; Wilton, 2001). We hypothesise that the under rep- 
resentation of dingoes from eastern and southern Australia in 
reference populations during both marker selection and analy- 
sis has confounded microsatellite-based dingo ancestry analy- 
ses. In light of these data, we urge wildlife managers to interpret 
microsatellite-based DNA testing data cautiously, particularly of 
dingoes from southern and eastern Australia. 

Genomic studies have identified historical or recent admix- 
ture in a wide range of species including bears (Cahill et al., 2015, 
2018; Kumar et al., 2017; Pongracz et al., 2017), hares (Jones 
et al., 2018), jackals (Galov et al., 2015), dolphins (Brown 
et al., 2014; Gridley et al., 2018), ducks (Stephens et al., 2020) 
and parrots (Qu et al., 2012). Hybridisation and genetic admixture 
have been an integral part of evolutionary adaption and specia- 
tion, but it also challenges wildlife facing urbanisation, changing 
climates and habitat fragmentation (Ottenburghs, 2021; vonHoldt 
& Aardema, 2020). Genotyping-by-sequencing or the use of tar- 
geted SNP panels may allow researchers to collect data from non- 
invasive samples like scat or hair for landscape scale monitoring 
of dingo and other wildlife populations. A key benefit of array or 
sequencing-based SNP genotyping is the ability to perform simul- 


taneous investigations regarding the genetics of biogeographic 
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variation, phenotypic traits, inbreeding, kinship relationships and 
demography, and facilitating more complete studies of dingoes 
and other wildlife. SNP-based genotyping with commercially avail- 
able microarrays offers the advantage of expandable data collec- 
tion, that is, SNP data generated at different times and locations, 
but with the same microarray, can be combined for a reproducible 
data set of high integrity. 

In a rapidly changing field, we highlight the challenge for sci- 
entists and wildlife managers to accurately identify admixture, di- 
agnose the processes driving hybridisation, assess the threat to a 
species or population and develop evidence-based conservation 
management plans to protect our wildlife. It is clear that as sci- 
ence advances, policymakers and wildlife managers need to rely 
on modern genome science and genome wide assessment methods 
to understand contemporary dog admixture in dingo populations. 
Furthermore, current policy and management plans regarding dingo 
conservation and management driven by microsatellite-based sur- 


veys should be revised and updated. 


4.4 | The identity and management of dingoes 


While dingo xdog hybridisation is not a pervasive threat to most 
dingo populations, we acknowledge that in parts of Australia, where 
there is a long history of suppressing dingo populations, levels of 
hybridisation are, by comparison, high. Thus, dingo x dog hybridisa- 
tion may be a threat to specific regional dingo populations. There is 
ongoing debate concerning when dog admixture in dingo popula- 
tions poses a threat to conservation and at what genetic threshold 
a dingo backcross becomes a dingo for purposes of conservation 
policy (Allen et al., 2017; Claridge et al., 2014; Crowther et al., 2020). 
We encourage policymakers to adopt a flexible framework such as 
that outlined by vonHoldt et al. (2018) which assesses local conser- 
vation goals, ecological function and a genomic synthesis of spe- 
cies concepts. At a minimum, it would be appropriate to expand the 
definition of dingoes in conservation policy and legislation to include 
historical dingo backcrosses using a genetic threshold of 93% dingo 
ancestry, as suggested by Allen et al. (2017). 

In Australia, the term ‘wild dog’ is widely used in policy and leg- 
islation based on the logic that the wild canine population is dom- 
inated by mixed ancestry dingoes or feral domestic dogs (Bird & 
Bowman, 2016; DEPI, 2013; Fleming et al., 2001; NWDAP, 2020). 
Given our demonstrated findings of limited domestic dog admix- 
ture in dingoes, future use of the terms ‘dingo’ and ‘feral domes- 
tic dog’ will more accurately reflect the identity of wild canids in 
Australia. Furthermore, such a shift in terminology aligns with calls 
from Australian Aboriginal people for government policy to better 
acknowledge and respect the value of dingoes as a native animal and 


culturally significant species (Costello et al., 2021). 
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